Analytic structure of the Landau gauge gluon propagator 
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The analytic structure of the non-perturbative gluon propagator contains information on the 
absence of gluons from the physical spectrum of the theory. We study this structure from numerical 
solutions in the complex momentum plane of the gluon and ghost Dyson-Schwinger equations in 
Landau gauge Yang-Mills theory. The resulting ghost and gluon propagators are analytic apart 
from a distinct cut structure on the real, timelike momentum axis. The propagator violates the 
Osterwalder-Schrader positivity condition, confirming the absence of gluons from the asymptotic 
spectrum of the theory. 



Introduction 

One of the fundamental properties of QCD is the 
absence of its elementary degrees of freedom, the 
quarks and gluons, from the physical spectrum of 
the theory. The associated problem of quark con- 
finement is a much debated issue [l| . In this discus- 
sion it is useful to distinguish between two notions 
of confinement. One is in terms of color confine- 
ment, i.e. the absence of colored states from the 
asymptotic, physical state space of the theory. The 
other is strictly related to the center symmetry of 
Yang-Mills theory. Both notions are not equiva- 
lent. If center symmetry is unbroken, there ex- 
ists a linear rising potential between static color 
charges in the fundamental representation of the 
gauge group, stretching out to arbitrary distances 
r of the charges. In QCD, the mere presence of 
fundamental dynamical charges breaks this sym- 
metry. Consequently, string breaking sets in at 
a sufficiently large separation R of the fundamen- 
tal test charges and the potential becomes flat for 
r > R. Thus, if confinement is defined in terms 
of unbroken center symmetry, QCD is not a con- 
fining theory Gluons are also not confined in 
this sense, since they life in the adjoint represen- 
tation of the gauge group. When two gluons are 
separated far enough from each other they pull two 
additional gluons out of the vacuum and dress up 
to form colorless bound states. This brings us back 
to the other notion of confinement: the absence of 
colored asymptotic states. A possible explanation 
of this absence are positivity violations. By defi- 
nition, the asymptotic Hilbert space T-L p h ys of col- 
orless physical particles has to be positive (semi-) 
definite, otherwise a probabilistic interpretation of 
its S-matrix elements would not be possible. Thus, 
positivity violations in the gluon propagator con- 
stitute a sufficient signal for the absence of gluons 
from the asymptotic spectrum of the theory. 

Another source of interest in the analytic struc- 
ture of the gluon propagator comes from Heavy- 



Ion collisions. Currently, there is great activity 
both from theory and experiment at RHIC and 
ALICE/LHC to shed light on the properties of the 
quark-gluon plasma (QGP), i.e. strongly interact- 
ing matter at large temperatures and/or density. 
Transport models like the Parton-Hadron-String 
Dynamics approach (PHSD) Q analyze the dy- 
namics of quarks and gluons in the QGP. An im- 
portant input into these calculations is the temper- 
ature dependent spectral function of the gluon, a 
quantity directly related to its analytical structure. 
Whether and how this structure changes below and 
above the deconfinement transition is currently an 
open question 

How can the zero temperature analytic structure 
of the gluon propagator look like? Based on studies 
of the gauge fixing problem, Gribov Q and later 
on Zwanziger [H[ suggested a form for the gluon 
propagator with complex conjugate poles at purely 
imaginary squared Euclidean momenta. A gener- 
alization with complex conjugate poles in the neg- 
ative half-plane of squared complex momenta has 
been suggested by Stingl in Ref. @ and has been 
recently explored in detail in the Refined Gribov- 
Zwanziger framework Ref. 0] ■ These forms have in 
common that they may pose a problem for the an- 
alytic continuation of the theory from Minkowski 
space to the Euclidean formulation used by lattice 
gauge theory and functional methods. Further- 
more, in their pure form they do not account for 
the perturbative, logarithmic running of the prop- 
agator in the large momentum region. An alterna- 
tive form with a branch cut structure for real and 
time-like squared momenta has been proposed in 
Ref. and found to compare well with numerical 
results for the propagator and its Schwinger func- 
tions in the Dyson- Schwinger (DSE) approach. 

All explicit calculations of the gluon propagator 
so far have been restricted to the real and space- 
like Euclidean momentum domain. Clearly, in or- 
der to pin down the analytic structure of the gluon 
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FIG. 1: Dyson-Schwinger equations for the gluon and 
ghost propagator. Filled circles denote dressed prop- 
agators and empty circles denote dressed vertex func- 
tions. 



propagator, an extension of these calculations into 
the complex momentum domain is highly desir- 
able. In this letter we report the first results of 
such a calculation. Within the continuum formula- 
tion of Landau gauge Yang-Mills theory we solved 
the coupled system of Dyson-Schwinger equations 
for the non-perturbative gluon and ghost propaga- 
tors in the complex momentum plane and extract 
the analytic structure at time-like momenta. As 
a main result we find analytic propagators every- 
where, apart from cuts on the real, timclikc mo- 
mentum axis. 

The framework 
The Dyson-Schwinger equations (DSEs) for the 
ghost and gluon propagators are shown in Fig. [TJ 
They form a coupled set of integral equations with 
renormalized, bare and dressed propagators and 
vertices. In Landau gauge, the explicit form of the 
propagators is given by 



D G (p) 



p^pA z(p 2 ) 



G( P 2 ) 



(1) 



with the gluon dressing function Z{p 2 ) and the 
ghost dressing function G(p 2 ). These functions can 
be numerically determined from their DSEs pro- 
vided explicit expressions for the dressed ghost- 
gluon, three-gluon and four-gluon vertices are 
known. Since these satisfy their own DSEs con- 
taining unknown higher n-point functions, in prac- 
tice one needs to truncate this tower to generate a 
closed and solvable system of equations. Certainly, 
meaningful results can only be achieved by careful 
control of the quality of such a truncation. 

A scheme which maintains multiplicative renor- 
malizability and transversality has been devised in 
Ref. for transverse projection of the gluon DSE 
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FIG. 2: Results for the gluon dressing function Z(p 2 ) 
from lattice calculations [12] compared to the result 
from DSEs 



this scheme is equivalent to the one used previ- 
ously in Ref. [k|. It uses a bare ghost-gluon ver- 
tex, a choice which is close to the results of corre- 
sponding lattice calculations [ll| , an ansatz for the 
dressed three-gluon vertex in terms of the propa- 
gator functions, see 0, [Io[ for details, and a van- 
ishing four-gluon interaction. Given this choice, 
the coupled system of DSEs can be solved numeri- 
cally. The resulting solution for the gluon dres sing 
function Z(p 2 ) has been discussed in Refs. [il. llCf 
and, for the convenience of the readers, is shown 
again in Fig. [2] together with corresponding lattice 
results (l2j . In the large momentum region, where 
the perturbative behavior sets in, both approaches 
agree very well. This is also true in the low mo- 
mentum region. In the deep infrared, the gluon 
dressing function displays the 'massive' behavior 
Z{p 2 ) ~ p 2 . Such 'decoupling' type of solutions, 



as opposed to 'scaling' [13|, have been suggested 



long ago 14 1 and have been revived in Refs. floT - 
17| . Large volume lattice results agree with this 
type of solutions [18[, although it remains a mat- 
ter of current debate whether problems with gauge 
fixing in the context of Gribov copies are already 
well under control lij- 22 1 . 

In the mid-momentum region around one GeV 
there are differences between the DSE and the 
lattice result on the twenty percent level which 
have to be attributed to the above discussed ver- 
tex truncations for the ghost-gluon, three-gluon 
and four-gluon vertex. Improvements for the 
drcssedghost-gluon vertex have been discussed in 



Refs. [23j, [24{. Furthermore, first studies of other 
types of ansaetze for the dressed three-gluon ver- 
tex are available [24[ and studies in the back- 
ground gauge Pinch-technique scheme emphasize 
the importance of poles in the longitudinal parts 
of the three-gluon vertex, see [25| for a review. 
While all these studies are interesting on system- 
atic grounds, the resulting solutions for the gluon 
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dressing function are not closer to the lattice result 
than the one shown in Fig. [2] The remaining differ- 
ence may therefore very well be attributed to the 
missing two-loop diagrams involving the four-gluon 
interaction. Indeed, pointwise agreement with the 
lattice data has been achieved within the frame- 
work of functional renormalisation group equations 
, where such contributions can be taken into ac- 
count due to the exact one-loop structure of the 
equations. Within the DSE framework the tech- 
nical complications arising from the two-loop di- 
agrams only allowed for phcnomenological treat- 
ments of these contributions so far [261 ] . and we 
therefore prefer to defer a study of the influence of 
these terms to future work. 

The numerical techniques necessary to solve a 
coupled set of DSEs in the complex plane have 
been explored up to now only in the context of the 
fermion propagator, see {27J and the appendix of 
[H . The basic idea is to shift the contour of radial 
integration in the loop integral into the complex 
plane such that singularities in the angular inte- 
gral are meliorated. For this work we adapted this 
method for the ghost and gluon system. Details 
will be given elsewhere. We have cross-checked 
our numerics also by employing a different method 
which solves the DSEs directly on a grid of com- 
plex momenta without any shifts in the integrals 
[28( . The results of both methods agree well for a 
large range of complex momenta. However, close 
to the timelike momentum axis, our first method 
clearly delivered much more stable results and is 
therefore to be preferred. 

Results and discussion 
Our results for the analytic structure of the gluon 
and ghost propagator in SU(N) Yang-Mills theory 
are shown in Figs. [3] and |U Let us first discuss 
the real parts of the propagators. In the lower 
diagram of Fig. [3] we see a spike structure of the 
ghost dressing function, which is located at the 
origin of the complex momentum plane. Since we 
have chosen a decoupling type of solution for the 
ghost-gluon system, the ghost dressing function is 
finite at p 2 = 0. At our numerical infrared cut- 
off |e 2 | = 10 -5 the value G(e 2 ) depends slightly 
on the direction from which zero is approached. 
In our calculation G(e 2 ) — 5 when p 2 — > + , but 
G(-e 2 ) = 5.005 T-i 0.004 when p 2 -)• 0" on the real 
axis. Thus, the real part of the ghost dressing func- 
tion is almost symmetric around the origin of the 
complex momentum plane. For the real part of the 
gluon propagator in the upper diagram of Fig. [3] 
the situation is entirely different. Again, the prop- 
agator is finite at p 2 = but shows large positive 
values for complex momenta and negative struc- 
tures close to the negative real momentum axis. 




FIG. 3: Results for the real part of the gluon propa- 
gator function D(p 2 ) and the ghost dressing function 
G(p 2 ) in the complex momentum plane including col- 
ored contour maps and lines. The displayed range of 
the gluon propagator is restricted in order to resolve 
smaller structures. See text for the extrema of ?RD(p 2 ). 




FIG. 4: Results for the imaginary part of the gluon 
propagator function D(p 2 ) and the imaginary part of 
the ghost dressing function G(p 2 ) in the complex mo- 
mentum plane including colored contour maps and 
lines. The displayed range of the gluon propagator 
is restricted in order to resolve smaller structures. See 
also Fig. [S] for the full scale. 

The positive spikes extend up to $tD(p 2 ) = 40 
GeV 2 staying definitely finite. Closer to the nega- 
tive real momentum axis the propagator becomes 
negative for \p 2 \ larger than some finite value on 
the negative real momentum axis. The correspond- 
ing narrow dip is finite in depth and sizably ex- 
tends 0.3 GeV 2 out into the complex momentum 
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FIG. 5: Results for the gluon spectral function and the 
ghost spectral function as a function of momentum. 



plane. The minimum value of the dip is approxi- 
mately 5i£>0 2 ) = -105 GeV 2 . 

Now let us discuss the imaginary part in Fig. [4] 
Here we clearly see a cut-structure emerging for 
both, the gluon propagator and the ghost dressing 
function along the negative real momentum axis. 
No further structure is seen in the complex momen- 
tum plane. We thus arrive at an important result 
of our study: the ghost and gluon propagators have 
nontrivial analytic structure only on the timclikc 
real momentum axis. This is in contrast to the ex- 
pectations from the studies of Gribov, Zwanziger 
and Stingl 0-0], which all assumed singularities 
away from the real momentum axis. We find no 
evidence for these. Within the present numerical 
accuracy, the cuts are sharply peaked but finite. 
Whether an even more precise treatment leads a 
singularity in ^sD(p 2 \, as assumed in the fits to the 
DSE results in Ref. [8j|, remains an open question. 

The cuts in the imaginary part of the ghost and 
gluon propagators are directly related with the cor- 
responding spectral functions, 

p G , g (p 2 ) = -Z{D G , g (p 2 )}/7r (2) 

with D G (p 2 ) = -G(p 2 )/p 2 , D g (p 2 ) = Z(p 2 )/p 2 
and the momentum p 2 on the negative real mo- 
mentum axis of the upper complex half plane. We 
therefore show them more closely in Fig. [5] Our 
numerical results are obtained on a grid of momen- 
tum points which arc displayed explicitly, whereas 
the interpolation is done via Chebychev polynomi- 
als and serves to guide the eye. For the ghost spec- 
tral function (and most of the gluon) this interpola- 
tion clearly works also on a quantitative level; how- 
ever, the interpolation for the gluon in the region 
0.5 < \p\ < 0.7 GeV can only be regarded as a qual- 
itative one. Whereas the ghost spectral function is 



dominated by the massless 1 jp 2 pole in the propa- 
gator, the gluon is clearly different. Rs imaginary 
part rises first, turns over and crosses through zero, 
turns again and then approaches zero from below 
at large timelike momenta. The exact locations 
and heights of the maxima in the positive and neg- 
ative region cannot be determined precisely within 
the present numerical accuracy. Nevertheless, the 
qualitative behavior is fixed from the explicitly cal- 
culated points shown in the plot. The gluon spec- 
tral function obeys the Oehme-Zimmermann nor- 
malization condition [29[ 

Z- 1 = J p g (s) ds, (3) 

where Z% denotes the gluon wave function renor- 
malization constant, with a deviation of 10 per- 
cent. This provides a measure on the accuracy 
achieved in the present computation. The nega- 
tive contributions to the gluon spectral function 
indicate its absence from the asymptotic spectrum 
of the theory. In general, the cuts in the ghost and 
gluon propagators signal the radiation of unphys- 
ical particles (ghosts and gluons) from unphysical 
particles. Moreover, the gluon is certainly not a 
massive particle in the usual sense. Nevertheless 
one may be tempted to define something like an 
"effective mass" m g for the gluon from the loca- 
tion of the positive peak in the spectral function. 
Within the present accuracy we find 

600 MeV < m g < 700 McV. (4) 

We stress again, however, that m g is not a measur- 
able quantity; strictly speaking, it is just the scale 
where positivity violations in the gluon set in. 

In this work we presented the first non- 
pcrturbative solution of the gluon and ghost prop- 
agators in the complex momentum plane together 
with an extraction of their respective spectral func- 
tions. We presented results for the decoupling case; 
a comprehensive comparison with scaling will be 
given elsewhere. Besides the considerable theoret- 
ical interest in these functions they are also a neces- 
sary input into the calculations of glueball masses 
within the framework of Bethe-Salpeter equations. 
Corresponding results will be detailed in a subse- 
quent work. 
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